Electronic origin of magnetic and orbital ordering in insulating LaMnOa 
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We derive a spin-orbital model for insulating LaMnOs which fulfills the SU{2) symmetry of 
S = 2 spins at ions. It includes the complete Cg and t2g superexchange which follows from 

a reahstic Mn'^"'" multiplet structure in cubic site symmetry, and the Jahn- Teller induced orbital 
interactions. We show that the magnetic ordering observed in LaMnOa is stabilized by a purely 
electronic mechanism due to the eg-superexchange alone, and provide for the first time a quantitative 
explanation of the observed transition temperature and the anisotropic exchange interactions. 
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The fascinating properties of doped manganites 
Ri-£cAa;Mn03, where R is a rare earth element, and A is 
a divalent element, were discovered almost half a century 
ago [Q, but the various phase transitions occurring un- 
der doping and in particular the phenomenon of 'colossal 
magnetoresistance' (CMR) are still not fully understood. 
The phase diagrams of Lai_a;(Ca,Sr)a;Mn03 show a 
complex interplay between magnetic, charge, and struc- 
tural order, so that all these ordering phenomena may 
affect CMR at least indirectly. It is therefore important 
to obtain first of all a full understanding of the mech- 
anism(s) stabilizing the observed order in the undoped 
insulating parent compound LaMnOs. This will be an 
essential element in putting together a satisfactory de- 
scription of the more complicated doped compounds, and 
recognizing which mcchanism(s) other than or in addition 
to double exchange |^ might be responsible for CMR . 

In this Letter we therefore reconsider the problem of 
the microscopic origin of the experimentally observed 
type of antiferromagnetic (AF) order in LaMnOa 
which consists of ferromagnetic (FM) planes ordered an- 
tiferromagnetically in the third direction (A-AF phase). 
As the magnetic order in LaMnOs couples to orbital or- 
der 1^] , one possible explanation might be the occurrence 
of a cooperative Jahn- Teller (JT) effect Q which induces 
a particular order of the singly occupied Cg orbitals 
However, while the JT effect plays a crucial role in charge 
transport |^ , we show here that a purely electronic mech- 
anism drives orbital and magnetic ordering in the man- 
ganites near the Mott- Hubbard transition [lC| ]. 

The local Coulomb interaction U is the dominating en- 
ergy scale in late transition metal oxides. If partly filled 
orbitals are degenerate, as in KCuFa or in LaMnOa, this 
leads to an effective low-energy Hamiltonian, where spin 
and orbital degrees of freedom are interrelated ]^,pd[| . In 
the simplest case of ions in KCuFa, such a model de- 
scribes spins 5 = 1/2 of Cg holes coupled to the discrete 



orbital variables. Finite Hund's rule exchange Jh re- 
moves the classical degeneracy of magnetically ordered 
phases pp^, and stabilizes the A-AF phase in conjunc- 
tion with the particular orbital order observed in KCuFs. 
Here we show that a similar state follows from a realistic 
S = 2 spin-orbital model for the c?^ ions in LaMnOa. We 
include also the t2g superexchange and the JT interac- 
tion and show that these, while unessential qualitatively, 
are important for a quantitative understanding. 

The superexchange between total spins S = 2 at the 
rf4 jyjnS-i- jQj^g originates in the large-JJ regime from vir- 
tual {cg or t2g) excitations, djdj # d^d^. A simplified 
approach proposed recently by Ishihara et al. [p^ empha- 
sizes the role of orbitals but violates the SU{2) spin sym- 
metry, and involves a Kondo coupling between Cg and t2g 
spins, which by itself is not a faithful approximation to 
the multiplet structure. The latter objection applies also 
to the model proposed by Shiina et al. In contrast, 
the spin-orbital model presented below follows from the 
full multiplet structure of the Mn ions in octahedral sym- 
metry, both in the d* (tL^g) configuration of the Mn'^+ 
ground state and in the d'^ and d^ virtually excited states. 

So we consider a spin-orbital model for the manganites. 



H — He + Ht + JT + Hr , 



(1) 



which includes superexchange terms due to Cg (He) and 
t2g (Ht) excitations, JT interaction (i/jx), and a low- 
symmetry crystal field (Hr). Our starting point is that 
each Mn'^^ (c?'*) ion is in the strong-field {t^gCg) Hund's 
rule ground state, i.e., the high-spin {S = 2) orbital dou- 
blet ^E. First, we analyze the strongest channel of su- 
perexchange, which originates in the hopping of an Cg 
electron from site i to its neighbor j. When we con- 
sider a bond oriented along the cubic c-axis, only a 
3^2 _ j,2 gigctron can hop and four d^ states may be 
reached: the high-spin ^Ai state [S — 5/2), and the lower- 
spin {S = 3/2) Ml, and M2 states (Fig. |lj). The 
rf^dj ^ '^i (^2g)^i(^2£(Sg) sxcitation energies require for 
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their description in principle all three Racah parameters, 
A, B and C e{^Ai) = A-8B, £(Mi) = A + 2B + 5C, 
eC^E) ~ A + 6B + 5C ei^) ^ A+UB + 7C. In view 
of the reahstic values oi B = 0.107 and C = 0.477 eV for 
Mn^+ (d^) ions one may use an approximate rela- 
tion C ~ 4_B, and write the excitation energies in terms 
of Coulomb, U = A + 2B + 5C, and Hund's exchange, 
Jh = 2B + C, parameters: e{^Ai) = U-5Jh, £(l4i) = U, 
eC^E) = U + ^Jh, £(^2) ^ U + fJn- Using the spin 
algebra (Clebsch-Gordon coefficients) and the reduction 
of product representations in cubic site symmetry ||l^] 
for the intermediate states, and making a rotation of the 
terms derived for a bond (ij) \\ c with respect to the 
cubic axes, one finds a compact expression. 
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where t is the hopping element along the c-axis, and ^^j^^ 
are projection operators for each bond (ij), 



'2P^CPK, (3) 



projecting on the orbital states, being either parallel to 
the bond direction on one site (P^.^ = 5 ~ ''"D 

: + t"), or parallel on 



pendicular on the other [Pj^ = ^ 
both sites. They are represented by the orbital operators 
r" associated with the three cubic axes (a — a, b, or c). 



q(6) 
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(4) 



where the cr's are Pauli matrices acting on the orbital 
pseudo-spins |a;) = (^o)' 1-^) ^ (i)' orbitals 



transform as \x) oc x — y and \z) cx (3z — r 



The spin operators Si in Eq. (g) are S — 2 spins, but 
otherwise resembles the spin-orbital model for ions 
in the cuprates . Both models contain superexchange- 
like couplings between spin and orbital degrees of free- 
dom. The orbital sector carries a discrete cubic symme- 
try, and is identical in both cases, while the spin problem 
fulfills the SU{2) symmetry, and different representations 
apply for the manganites {S — 2) and for the cuprates 
(5 = 1/2). We emphasize that the Hamiltonian He is 
not equivalent to that of Ref . jl^ in any nontrivial limit. 
A common feature is that FM interactions arc enhanced 
due to the lowest excited ^Ai state, but the dependence 
of the magnetic interactions on Jh is quite different, and 
it gives a different answer concerning the stability of the 
A-AF phase. The balance between AF and FM interac- 



tions is also different from that in Ref. |14| due to the 
multiplet structure of Mn^+ . 



A similar derivation gives the t2g superexchange [iq ], 

where Jt = (Jn + J22 + J12 + >/2i)/4. The exchange 
elements, J,„„ = tj/£(^r,„,4r„), where = t/3 is 
the hopping between the t2g orbitals, result from local 
dfd'j ^ c^f (^2g 69)^^(^29 6g) excitations within a (ij) bond, 
with energies ei^Ti,'^Ti) U + SJu/i, e{'^Ti,^T2) ~ 
U + 2Jh/S, e{^T2,^Ti) ~ U + AJh, £(^^2,^1^2) ^ U + 2Jh, 
where ^J"^ (^TVt) stands for the symmetry of c?f (d^) ex- 
cited configurations, respectively. 

The manganite model (Q) is completed by the JT term 
which leads to static distortions and mixes Cg orbitals , 



H 



JT 



27'f.^ 
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with 'P^jj-^ — 2Pi^Pj^, and by the tetragonal crystal field. 



(7) 



The strength of Cg and tg superexchange can be esti- 
mated fairly accurately from the basic electronic param- 
eters for the Mn ion as determined from spectroscopy 
|]l7| , ^ with an estimated accuracy of ^ 10%. We thus 
use U = 7.3 eV and Jh = 0.69 eV, and taking into ac- 
count that the Mn-Mn hopping occurs via the bridging 
oxygen, t — 0.41 eV as follows from t = t^^/ IS. with 
Mn-0 hopping tpd = 1.5 eV and charge transfer energy 
A = 5.5 eV @. This yields J = t^ /U = 23 meV and 
Jt — 2.1 meV. The accuracy of these parameters may 
be appreciated from the resulting prediction for the Neel 
temperature of CaMnOs, where a similar derivation gives 
Ht - 2Jt{^Si ■ Sj - 1) in terms of Mn^+ spins S = 3/2 
and Jt ~ Jt(l + Jh/U). With our present estimates we 
obtain Jt = 2.3 meV and thus Tjv = 124 K, in excellent 
agreement with the experimental value Tj\j = 110 K 

When considering the manganite (d'*) model (^, it is 
instructive to treat Jh and Ez as freely variable param- 
eters in order to appreciate the physical consequences of 
Hund's rule multiplet splitting and orbital degeneracy. 
The cuprate (d^) model exhibits symmetry breaking into 
classical states with simultaneous spin and orbital order 
|H,|2|, and similar behavior is expected here |2|]. We 
have considered classical phases with two and four sub- 
lattices, and mixed orbitals (MO), \ifia) — cos6i\ixa) + 
sm6i\iz(j). The mean-field (MF) phase diagram of the 
Gg-part of model (|l|), = + Ht, at T = is similar 
to that of the cuprate spin-orbital model at large 

positive (negative) E^, one finds AF phases with either 
I a;) (AFxx) or \z) (AFzz) orbitals occupied, while MO 
phases with orbitals alternating between the sublattices 
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{6i = ±9, with cos 20 < 0) are favored by increasing Jh- 
If Ez < the spin order is FM (AF) in the (a, b) planes 
(along the c-axis) in the MOffa phase, while at i?^ > 
two similar phases, MOAFF and MOfaF, are degenerate. 
For the parameters appropriate for LaMnOa {Jh/U ~ 
0.095) one finds a MOffa/MOAFF ground state, i.e., A- 
AF magnetic order, while a FM (MOfff) phase is found 
only at Jh/U > 0.12. The region of stability of the A-AF 
phase is modified by t2g-superexchange [Fig. ||(a)], but 
this change is small as Jt ^ J. Thus, the observed A-AF 
magnetic order in LaMnOs is caused by the orbital de- 
pendence of the Cg-superexchange and not by competition 
between FM Cg- and AF <2£rsuperexchange as proposed 
in Ref. [|3| (where an unrealistically large Jt was used), 
supporting the qualitative results of Ref. [0 . 

Although the MF phase diagrams are modified signif- 
icantly by JT coupHng [Fig. ||(b)], the A-AF phase sur- 
vives around Jh/U = 0.095. In fact, the JT interaction 
(^ by itself enforces alternating orbitals with cos26' = 0, 
which favors AF spin order, thus stabilizing at small 
Jh/U the MO AAA phase, promoted further by finite Jt- 
But at larger Jh/U, even though the JT interaction sets 
the stage by inducing orbital order as such, the actual 
magnetic (A-AF) and orbital (cos 26* ^ 0) order are en- 
tirely due to the Cg-superexchange interactions (^. 

Finite temperature behavior was investigated in MF 
approximation, with (a), (S), and (aS) constituting in- 
dependent order parameters |2|]. As the largest interac- 
tion is in the pure orbital (a) channel, one may estimate 
the JT coupling k from the temperature of the structural 
transition, T^^p w 750 K. The electronic interactions 
contribute ~ 440 K (Fig. |), and the rest, 6k ~ 760 
K, comes from the JT term [Q. Thus k ~ 11 meV, 
and we have adopted the representative value k/ J — 0.5. 
We then calculated the temperatures Tc for the possible 
magnetic transitions (Fig. ^), taking into account that 
orbital order with cos 26* ~ already exists below Ts, 
and calculating selfconsistently the corresponding order 
parameter (a) at finite T. The spin order sets in simul- 
taneously with a modification of orbital ordering towards 
cos 20 7^ 0. We find that the preexisting structural tran- 
sition plays an important role at finite T and reduces the 
magnetic transition temperature, being otherwise Tc J 
. The results are consistent with the phase diagrams 
at T = (Fig. H), as the magnetic transition corresponds 
to the same order as found at T = 0. For the A-AF 
(MOffa) phase we find ~ 106 K [|3|, in reasonable 
agreement with the experimental value of 136 K p4[ . 

The magnetic interactions in the A-AF (MOffa) phase 
may be found using averages (T'".^^) of the orbital projec- 
tion operators (^) at = 0. They are FM in the [a, b) 
planes iJ(a.b))i a-nd AF in the c-direction (J^,) (Fig. |^). 
Both large Jh/U and k — 0.5 J play a decisive role in de- 
termining the actual composition of the orbitals, and we 
find J(a,b) = —1.15 and Jc = 0.88 meV, somewhat higher 



than the experimental -0.83 and 0.58 meV [^^. However, 
their ratio, Jc/\J(a,b)\ = 0.77, agrees very well with the 
experimental value of 0.7 — 0.72 [^^; in contrast, it would 
amount to 1.04 for k = 0, and to 2.25 if in addition the 
orbitals were chosen to satisfy cos 20 = —0.5. 

Summarizing, a coherent overall picture has been ob- 
tained for LaMnOa, which is even quantitatively satisfac- 
tory. It includes simultaneously the full (e^ and t2g) su- 
perexchange and the JT effect, and shows that the orbital 
dependence of the e^-superexchange, a purely electronic 
mechanism, is responsible for the observed A-AF order. 
We emphasize that no fitting of parameters was needed, 
and the used values of Jh, U, and t, known with an accu- 
racy of ~ 10%, allowed to deduce the value of the JT cou- 
pling K, and gave Tc, J(a,b) E^nd Jc within 30% from the 
experimental values. We thus believe that Hamiltonian 
(|l|) provides a realistic starting point for understanding 
how the delicate balance of magnetic and orbital inter- 
actions in LaMnOa is affected by doping, leading to a 
change of magnetic order and to the possible onset of an 
orbital liquid state j2^ . 
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FIG. 1. Virtual d^d^ excitations which generate 

effective interactions for a bond (ij) \\ c-axis: (a) for one \x} 
and one \z) electron, and (b) for two \z) electrons. 

FIG. 2. Classical phase diagram of the manganite model 
(|): (a) no JT effect (k = 0), Jt = (full lines) and 
Jt = 0.092 J (dashed lines), with the AFxx and AFzz phases 
separated by a MOAAA phase; (b) including JT effect 
(k = 0.5 J), Jt = (dashed lines) and Jt = 0.092 J (full lines). 

FIG. 3. Magnetic transition temperatures Tc (J = 23 meV, 
=0, Jt = 0.092 J, K = 0.5 J) for: MOAAA (dashed fine), 
MOFFF (long-dashed line), and MOFFA (full line) phases, 
and Ts for the structural (MO) phase transition at k = 
(squares). The dotted line indicates realistic Jh/U = 0.095. 

FIG. 4. Exchange interactions Jab and Jc in the ground 
state for increasing Jh/U, for J = 23 meV and: Jt = 0, 
K = (dashed lines), Jt = 0.092 J, n = Q (long-dashed lines), 
and Jt = 0.092 J, k = 0.5 J (full lines). The inset shows cos 26. 
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